UNCLASSIFIED 


AD NUMBER 
NEW LIMITATION CHANGE 
TO 


Approved for public release, distribution 
unlimited 


FROM 
Distribution authorized to U.S. Gov’t. 
agencies only;. Other requests shall be 
referred to ASD, Wright-Patterson AFB OH. 


AUTHORITY 


USAF ltr 4 May 1977 


THIS PAGE IS UNCLASSIFIED 


ASD/XR-TR-75-1 


ASDIR-1TI 
VOLUME 11 
PROGRAM DESCRIPTION 


-- 


January 1975-7 


~ 
a! 
i <> _~ Distribution Limited to US Government agencies only; Aeport 
' SOD  eontains test and evaluation information; January 1975. 
i Other nequests for this document must be neferred to ASD/XRH, 
Ly Wright-Patterson AFB, OH 45433. 
ae 


es 
an) Copy available to DDC does not 
- yommit fully legible reproduction 


‘ DEPUTY FOR DEVELOPMENT PLANNING 
AERONAUTICAL SYSTEMS DIVISION 
WRIGHT-PATTERSON AIR FORCE BASE, OHIO 45433 


NOTICES f 


When Government drawings, specifications, or other data are used for any 
purpose other than in connection with a definitely re1ated Government pro- 
cureiment operation, the United States Government thereby incurs no respon- 
sibility nor any obligation whatsoever; and the fact that the Government 
may have formulated, furnished, or in any way supplied the said drawings, 
specifications, or other data is not to be regarded by implication or 
otherwise as in any manner licensing the holder or any other person or 
corporation, or conveying any rights or permission to manufacture, use, or 
sell any patented invention that may in any way be related thereto. 


This technical report has been reviewed and is approved for publication. 


Colonel, USAF 
Deputy for Development Planning 


UNCLASSIFIED 
SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) ‘ a tt 
READ INSTRUCTIONS ; 
REPORT DOCUMENTATION PAGE pee ee aS ES te rage °° 
Ji. REPORT NUMBER... meme siog hia tes 2. GOVT. ACCESSION NO. 3. RECIPIENT'S CATALOG NUMBER 
ASD PXR=75-1 —Voligeeti— 
- [4. TITLE Cand Subtitie) : DERREORT BPE RIOD COVERED 
SO ASUIRET Sete ‘of CS Mad oese hae 
\* || Volume IL, \ July 72 ~ZDec 74, | 
PROGRAM DESCRIPTION , PERFORMING ORG. REPORT NUMBER 
. = — 7 8. CONTRACT OR GRANT NUMBERS) 
(po \ Charles W./Stone/ Ist Lt., USAF PIG! 
“~.f} Stanley E. jTate f . \/ td Ol = Z 
} Rene nag: N/A Per fa bsos, 
9. PERFORMING ORGANIZATION NAME AND ADDRESS 10. PROGRAM ELEMENT, PROJECT, TASK 
peput for Deve lopment Plarming AREA & WORK UNIT NUMBERS 
ASD ARH s 


Air Force Systems Command - USAF 
Wright-Patterson AFB, Caio 45433 


11. CONTROLLING OFFICE NAME AND ADORESS wo | 42. REPORT OATE 


De for Development Planning Tanger 75 _— 
ASB AR © 3 stems Co eet |" LA PIS NUMBER OF PAGES ————J 
ae ght-Patterson Apemape. Aree 293 | 


15. SECURITY CLASS. (of this report) 


UNCLASSIFIED | 
DECLASSIFICATION/ DOWNGRADING 


iSe. FH 
| SCHEOULE 


N/A 


4. MONITORING AGENCY NAME & AOORESS/(if different from Controlling Office) 


36. DISTRIGUTION STATEMENT (of thia Report) 


Distribution limited to US Government agencies only; report contains test and 
evaiuation information; January 1975. Cther requests for this document must 
be referred to ASD/XRH, Wright-Pattérson AFB, CH 45433, 


ets 


17. DISTRIBUTION STATEMENT (of the abstract entered in Block 20, If different from Report) 


H 18. SUPPLEMENTARY NOTES 


i 


19. KEY WORDS (Continue on reverse aide if necessary and identify by block number) 


Infrared Signature Infrared Detectability 

Infrared Signature Prediction Infrared Suppression 

Computer Modelling Infrared Analysis | 

Infrared Survivability ASDIR 
¥ 


20. ABSTRACT (Continue on reverse side if necessary and identiiy by block number) 


The Aeronautical Systems Division Infrared Signature Prediction Program 
(ASDIR) is a state-of-the-art functionally modularized computer code. 
Volume #{ contains a fimctional and technical description of each of the 
programs used in ASDIR as well as a complete program listing. 


DD en, 1473 E0ITION oF 1 NOV6S IS OBSOLETE UNCLASSIFIED 
ae a _ SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 


— a oft : oa ie 
Ve 5 fr fot ee" wee 


TABLE OF CONTENTS 


GENERAL PROGRAM DESCRIPTION . 2. ee ee ee eee ee tee we es 


PROGRAM CONTROL 
DSDIRG 2. oe ao ee ee 6 BE na, BE) we et BE ws eter A ene 


Cr 2 


STGNIR: sorte sel fa) we Sole eran ere cel er GE, Wt Way ate a, OG SE le 


REFERENCES 


Ss 6 © © © © = © © © © # © © © © © e © © © 8 2B es 


APPENDIX - PROGRAM LISTINGS .. 2... 1... a eee Baa ie sett ABE 


OVERLAY (XRHP, 0, 0) 
DEDURD. cin, otal tee Se et i we, ah Se Hoke Gee es ae ed 


eS ee a eee es ee a ee ait 


i ae a el ee a 


HUNG Nay a8 eee bites Spe ete Monies. ve 2s holds Gag eee, aries Set he 
MRIS ses Sead. abate See Sates ak ei igh gp Sec Ne a Ges eM Se et aes eee SASS 
GRUER 8. Be. cd ote co: eh) Ge ye duet vic Gerais ey tal bese 
COND OF dine GING Bes Us RASA ha wel oe go ia Ne BIE Be gece Ce aL Whe D 
TOA EIN yt ete ge cpt boo agentes, Gal poy abs) Ja6 al det bee gee Ge Seiden oe? teste God per eG ate as JAR 


OVERLAY CXFHR, 2, 0) 
TW shave ge ag Bs Gh eS og «a ogg ob ea Be Swans Soe a a 


OVERLAY (XRHP, 2, 0) 
PACU cor geod Vacate hn de ue tae hhh Se wsde ok nee AED A- 
VON olor 2 Gene val’ dh co cee ah eee Me AE gh ocak Ca Re len tas i 4 A 
PIURAS Ale ak ce oe Sivig gil eed eae a nde aia a, tas ow 
INTERD A 


ET ae. od ae, Se Gy ta. Gone ey ae a, eye OS Bt a Le A decide dh, Yet eek doje se . A-63 
AVE... 
SETTAU 
“LAUCAL 
PEANGKS 2°; ecsy eras ah te sath aia Stee “ates aoeol. gah Sadar Bt Satay fest fe » + A-75 
TNTERPO® 002: os: eer an Sateen Sater Sad es SS - 
TAUNZO «6 6 se ww ew we ee te tee 
ERE stay Mot 2 Mee! baa et eee ee aay? 


OVERLAY (XRHP, 3, 6) 
DATENE <6: ha) store? Bt ces 2 és ‘ay To we 


OVERLAY (XRHP, 4, 0) 
PIMNEE sou eg get er ahaa See: Beet Me fan See fe ode geek eon Od 
BUN so ate Ge ete ete ee eens ee ins ete ta, gies ti ABE 
BEIM wise) 50 deve, Gees ee Racal we hw ca uk he ke eh odin ao ah, SACO 
MIG of checig cates Sah alee Sh tated hale © Rae aden eet AA ie views » AS 
OVERLAY (XRHP, 5, 9) 
DERM a s06 Se GP Sine wane ob eubie de atid is ee eer ae A-98 
. DIMPENT 24 0 dedi: 00 Sod. done beta boot Sr ave me pe ceccatne a) WO 
i PIIPLE “isasxte-a vk On ek ore ' A-101 


ok a dh ae ee en es 


OVERLAY (XRHP, 6, 0) 
BEER: See eet theres: Se fa BoB ee et cone tals wanna ont: es HOG 


OVERLAY (XRHP, 7, 0) 
MELO tases of base ta ads feo Mau ade eo ee re ree ae ss Mee 


OVERLAY (XRHP, 11, 0) dl 
REMUN <b be de Wy ada Beam eds Ge Ge hs oe ahs eR ee, AO ; 


| 
{ Ld 
mil 
4 
RING aoe ots Toes tes Montes aes eh howto: ca tes as He Se anes, Bt eee ld A-120 
PRAT Riese Ser haha EUG wats Se ep rete hae eh, Gt Ege ets ee ae Asy2t 
OREM icc es del cee de Scot eg hoe Aas sn ge ce ma oe els A-122 
STORE caxns< enced ol, be aw end ear ates te ES gin Mo eg Far A-124 
STONE. ver ges tech Boa amide, Bogie et edt ae ae fw Pi an nk A-125 
RARE DS 5p ts oe woe Bh se ee Gee Bre lel eene oe te ie A ee A-126 
GRIER ol o.ta: 4 Sat es ee ge Se ae ae wy Adoni eee 2 A-127 
PAG coo eo At oo ak ete cae he asthe) eda the, eae eed A-128 
PRNIE e saretee tan cya Ve ahaa teeta oe “bn Tee as tah een anh . A-129 
OVERLAY (XRHP, 11, 1) 
PANES opethy oo ok toe ie es Geen orate eta eee tee ee ne Ee ae A-130 
OTRO) diae a RD  iig tercel cine 0, fe ye Wiad eat aie tate ete Aad a, A-133 
OVERLAY (XRHP, 11, 2) 
2 1, OR a Ree ae ht oo Re ea re A-134 
ETE ca hir pakorse eicct hue? AP eee Glas te et ee at: Ab ee we ah, etapa tage oe A-138 
QVERLAY (XRHP, ll, 3) 
CONE 9 ou as et po dy, Oo ete Ae ed Ba ett en Ses A-139 
DROUIN enka co Pk cs eae a aS a eo ee De ne A-147 
OVERLAY (XRHP, 11, 4) 
CALFOR . oe @ e re «© @© ©» 8 © 2» © 0@ © © @ @ © © @ 9 . ee © @ © A-148 : } 
OVERLAY (XRHP, 11, 5) 4 
HOIPT 4 oe se Sie Bea de Ge Pe recccon ae pele eta me REISS qi 
AVERHT . 8 @ ee . . . ee © © 8 © © @ © * © © «© © © © © # #© @ @ A-151 f | 
HEAT ess tat Mao ae ch xt ee. RR ah ae Ba rer ee ee rene > bea 6 AEVSS iz 
PUIRBLTS o sr ccs an prvees tee erase By Tah as eee cay ee ee cae spend. ot wee AW154 {4 
BEATER hd vk Actes Bioae tgs. cieeli Sk webu Solar ah Te th ooh os tates x A-159 | 
VERE S Me ils een h ce, 409e. Seren ce Ceoaahies & ati ao ewig hsb evoke ete A-160 Vi 
OVERLAY (XRHP, 12, 0) 
VEEW Se cena ng. Se tape cis We hie ah EVE <x tA sete ae Fame ce calebh ME % A-162 ea 
PREC a9 135 econ Ca irda th re ee a wae A-176 , 
UND oe andi o> as eae ay aad. a, Ile: Aran Ge Serva ate wee tanta A oe ps A-177 4 
COB IG Wilh rh Bes rane eee ee ara tance aan ey eee we ad A-178 : 
MTS a eigrey Gur gue as Rtee ay aa ae eh se ete A-180 
MES agit a8 anh Oe ia OTe 6 he, as ee eee) Une ; A-181 } 
BONDE iceeeke, cerca seta ee ee A bes cesar atari an aa Ot scent A-182 
CODE ona rks cy Ste aa Sah bt a tke ed, Gace beaten! A-183 
CANBEES Scan ta Sa eet as ce clang Br alge sa Ser yuh el-vien aage lees A-184 
BOE ep gic eae Ss th pecs nee ans Soa ae este cae aetaae ee carer: eA: Mae A-186 
OHI, he oh 88 de Ge Ant he nd oo i es ies che ncaee HE Binh se eg A-188 
RUNGE sat ses este hash Oar tas Bac en oe ae ea Seer Se ea a card Sea hed A-189 
BUCHAN oda laa) 60S 9h Wie eer a oi Ragens ade an nama A-190 
PU CMUN RR eerie nage Bho elect ir cee taba, gh ihn gas ata a ae A-191 
EB SING: Ae yghtasiactas cause te oooe ve te! Cte Bea Wena & A-192 
OVERLAY (XRHP, 13, 0) 
: PAST? eck ween eat wate Stak ee We cee. ee aa EA Gh Oe a eek a A-193 ; 
cs CRONE fy se arraniceet Ante BGG cs ag Atk odes eae Sctls ety eel ok ied ate A-196 : 
MAL POR score Sage ety ee ee: Sik es nae, hee eo ote ae eae tak eg A-197 ' 
MAT GYM co et ce Soke ea" goo eS Able kas, Sane ene de ates Boek A-198 { 
i i 
I 
: 
: ‘ ' 
a 


DINVRT 
WRIT14 
KEDIT . 
CEDIT . 
SCRPTF 


OVERLAY (XRUP, 


TAIL. 
PLOTT . 


, 


OVERLAY (RHP, 


TATLT . 


RITE. 
VECANG 


OVERLAY (XRHP 
TAIL? 


» es © © es © 


14, 


14, 


. 


0) 


A-199 
An 201 
A-202 
A-205 


Ax 204 


A-205 
A-290 


A-298 
Ar 209 
A-210 
A-245 
A=216 
A-217 


A-218 
A-219 
A-22 

A- 222 
A-224 
A- 225 


eee 


ot 


ASDIR-II PROGRAM DESCRIPTION 
EE EE EOCRIPTION 


INTRODUCTION 


This Program Description outlines the structure of ASDIR-II in order 


to give the users the essential knowledge required for program implementa- 
tion, use and modification. Simplified structural and flow diagrams of 
the major programs and short functional descriptions of ail programs and 
Subroutines are given to help describe ASDIR. 


The level of detail presented 
in this manual is consistent with the availability of reference and docu- 
mentation material, program difficulty or progr i 
tional accuracy, 


Essentially this manual is wr 
Sequences in ASDIR, with modificat 
and clarity dictate. 


itten in the same order as the calling 
ion only when structural considerations 


Please note, this manual is intended as an aid to the user to allow 


better understanding of each program's function and not as a users Manual 
(see ASDIR-IT, Volune I - USERS MANUAL). 
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GENERAL PROGRAM DESCRIPTION 


The Aeronautical Systems Division's InfraRed Signature Prediction 

Model (ASDIR) is a computer code utilizing the overlay structure shown 
in figure 1. These overlays are organized into three major functionally 
modularized computational groupings, figure 2: (1) Hot Parts Analysis; 
(2) Plume Gas Dynamics; and (3) Signature Prediction. A fourth overlay 
element establishes the required data initialization and input/output 
operational capabilities. 


This volume is organized in the following functional format: 


. ASDIR Z - Pro::xan Control 

Data Initialization § I/O 

. Hot Parts Analysis 

Plume Gas Dynamics 

. Signature Prediction 

Program Listing in Overlay Form 


Aum PWN Re 


PROGRAM: 


FUNCTION: 


INPUT: 


OUTPUT: 


SUBROUTINES : 


DESCRIPTION: 


ASDIR 2 (Page A-4) 


This program is the main controlling routine of 
ASDIR (the Aerenautical Systems Division InfraRed 
Signature Prediction model). ASDIR's prime 
function is to provide a logical sequence of calls 
to the functional programs. 


None 
None 
All overlays and data blocks 


See figures 1 and 2 and program listing on Hi 
page A-4, ; 
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PROGRAM: 


FUNCTION: 


TNOUT: 


SUBROUTINES: 


DESCRIPTICN: 


INPUT (Page A-43) 


Provide a vehicle by which the desired calculation 
parameters, cose definition and data initialization 
can be read into ASDIP. 


7 


All contro] ané comrutational variaodles in the 
CASE nameiist. 


Computationaily compatible case and dimensional 
data. 


None 


for complete 
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PROGRAM: 


FUNCTION: 


. INPUT: 


OUTPUT: 


SUBROUTINES : 


DESCRIPTION: 


DATINT (Page A-79) 

To initialize the calculation parameters. 
None 

Initialized calculation parameters 


None 


This program has been written to initialize the 
common block data to the default case. See 
ASDIR-II, Volume I, User's Guide for default 
options. A secondary function of this program 
is to allow use of ASDIR on systems in which a 
zero core initialization is not available. 
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PROGRAM: FILTER (Page A-106) 
FUNCTION: Provides a means of accessing stored filter 
' response dat: or ontering a filter response 
of the user's choice. 
TNPUT: Filter number or (ilter hand width and ~esponse. 
OUTPUT: Filter response. 
SUBROUTINES: he 
DESCRIPTION: Stored filter responses are in data statements 
with modific.ition possible through namelist FILT. 
: Tais program also can ve used to initialize the 
2 wavenumber culculation intervals. 
‘ 
{ 
% 5 
3 
. q 


PROGRAM: 


FUNCTION: 


INPUT: 


OUTPUT: 


SUBROUTINES: 


DESCRIPTION: 


XPLOT (Page A-109) 


To provide a means of file manipulation for 


either punching or printing spatial signature 
data. 


Calculated spatial signature data. 
Punched or printed spatial data. 


None 


This program requires the use of a punch file 
and an addition file TAPE 8. 


PROGRAM: 


FUNCTION: 


INPUT: 


AT PIE: 


SUBROUTINES : 


DESCRIPTION: 
PROGRAM: 
FUNCTION: 
INPUT: 
OUTPUT : 
SUBROUTINES: 
DFSCRIPTION: 


Block Data H20 (Page 


Data Rlock for water 


None 


A-15) 


vapor band parameters. 


See References [1, 7, 3} 


Block Data COQ2 (Page A-7) 


Data Block for carbon dioxide band parameters, 


None 


None 


None 


See References [1], 2, 
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PROGRAM: 


FUNCTION: 


INPUT: 


OUTPUT: 


SUBROUTINES : 


DESCRIPTION: 


SIGNIR 


To calculate the hot parts signature and 
descriptive parameters required in ASDIR. 


Engine description, gas properties and case 
definition (See ASDIR-II, Volume I, User's Manual). 


Primarily the information which is transferred from 
the SIGNIR link to ASDIR is an equivalent black body 
area and equivalent biack body temperature. 


See Figure 3. 


With the exception of the modifications required 

to adapt SIGNIR to the ASDIR format the hot parts 
calculations are performed with the original LIV/ 
SIGNIR-I computer code. This code was chosen 

for only one reason; it has proven to be the best 
available methodology to predict engine hot metals 
Signatures regardless of program size, operating cost 
or complexity. 


The description of the infrared signature prediction 
program (SIGNIR) is provided by the following 
discussions which cover a general description of the 
program, the program computational procedure and the 
analytical methods utilized in obtaining the solution. 


The program SIGNIR is a digital computer program written 
for the purpose of providing predictions of the hot 
metal infrared emission from aircraft engine exhaust 
systems. The program is designed to be applicable to 
axi-symmetric turbojet, turbofan, or turboshaft engine 
exhaust systems. It predicts the spectral intensity 

of the radiant energy emitted from exhaust system hot 
metal in the wavelength band of from one to 15 microns. 
In general, the information required by the program is 
as follows: 


. Exhaust system physical characteristics 
. Engine operating conditions 
Special surface cooling flow conditions 
Exhaust system surface properties 


The predictions provided by the prog: un for the 
combination of a selected maximum of 20 engine off-axis 
angles are: 
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Spectral radiant intensity nonattenuated by the atmosphere 
Selected radiant intensity bandwidth summations nonattenunted 
by the atmosphere 


Optional exhaust system information which can be requested from the 
program is: 


. Internal fluid flow properties 
. Surface boundary layer data 

Internal and external geometric radiation view factors 
- Surface temperature distribution 
. Thrust information 


The process of computing the intensity of the infrared energy emitted 

from aircraft exha:.st system hot metal requires a knowledge of the 
metallic surfaces temperaturc, emissive properties and their relative 
physical positioning within the system. From this information, the radiant 
energy emitted and reflected from the éxhaust system to a point in spiice 
can be determined. The canputational procedure established by this program 
for the prediction of exhaust system hot metal infrared radiation is shown 
by major subjects in the general block diagram presented on figure 3. 


The input data required by the program provides a physical model (nodal 
system) of the exhaust system configuration, the properties and conditions 
of th: exhaust gases entering the system, band width and location relative 
to the exhaust system. The program undertakes a one-dimensional, isentropic, 
Compressible flow analysis which generates the axial distribution of exhaust 
gas flow conditions and properties within the system. Where two gas 

streams at different temperatures flow within a single passage, compound 
flow considerations are utilized. The program adjusts the flow ut the 
nozzle throat to balance with the ambient pressure or a chocked condition 
dependent on the system's entrance flow data. 


With the internal flow properties defined, a surface boundary layer analysis 
is conducted to obtain the coefficients of convection and friction along 
each surfuce within the system. This analysis is primarily dependent on 
the surface's axial static pressure distribution and the upstream boundary 
layer conditions. An area weighted average coefficient of convection for 
each of the model's surface nodes (isothermal regions established by the 
input data) is computed for the system thermal analysis. The uxial force on 
each surface, intemal to the system, is computed utilizing both the 
frictional and pressure-area forces. 


Descriptions of the routines associated with SIGNIR are presented in Table 
1. Further definitions of the routines can be found in ASDIR-II, Volume 
III, Reference Documentation or in Reference [4]. 


The program STGNIR, other than a few minor programming changes, remains 
as developed by LTV. The only link between SIGNIR and ASDIR-II is Common 
Block OUI, containing the following information: 


iquivalent Black Body Area, cme 

Equivalent Black Body Temperature, deg-K 
. Viewing Aspect Ancles, degrees off tail 
. tlumber of Aspect Angles 


ect eo er Se 


ROUTINE 


NAME 


ARCCOS 


ARRAY 


AVER 


AVERFS 


AVERHT 


AXAREA 


BODY 


CALFOR 


CEDIT 


CHECK 


TABLE 1 


DESCRIPTION OF SIGNIF ROUTINES 


BASIC FUNCTION 


Determines the angle, in radians, correspond- 
ing to a given cosine valuc. The angie will 
be from tero to If/2 radians. 


Describes in two-dimensional coordinates the 
inner and outer surfaces of the flew streans 
using the surface node coordinates. 


Averages the cooling wall temperatures for 
the surface nodes along the cooicd portion 
of the configuration surfaces; prints these 
averages, 


Determines the average values of the parameters 


necessary for force factor (pressure area 
and friction forces} calculations for each 
surface node along the configuration surfaces. 


Determines and prints the average heat transfer 


coefficient for each surface node along the 
configuration surfaces. 


Determines tne cross-sectional area of the 
fluid flow streams as a function of the 
configuration axial distance. 


Determines for the internal view factors, if 
the iine of sight between points on two nodes 
is tangent to a shadowing body. 


Prints the average force factor for each 
surface node and calculates and prints the 
total net surface force factor for the 
configuration. 


Checks, for internal view factors, the visi- 
bility between two nodes to determine if 
shadowing may occur. 


Selects the convection and conduction heat 
paths from input data and forms the conduction 
matrix. 


Checks the data returned from routine BODY 
for extraneous solutions. 
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PAGE 


A-176 


A-38 


A-118 


A-151 


A-126 


A-186 


A-184 


A-203 : 


A-188 


COFLOW 


COND 


CONFLM 


CPAGMA 


CRODEL 


CUBIC 


DINVRT 


DIS 


FEDIT 


FILMCL 


TABLE 1 (Cont) 


BASIC FUNCTION 


Analyzes the results provided by routine 
FINDIT to describe the nature of shadowing 
between two nodes. 


Calculates the compressible flow properties 
for the flow streams as a function of 
pressure. Adjusts the flow for cooling or 
ambient pressure at the exit of the 
configuration. 


Reads in those values of heat transfer and 
fluid lump temperatures for nodes whose 
thermal paths are not connected with the 
fluid streams. 


Calculutes a pressure balance and surface 
temperatures for a surface cooled by convec- 
tion-film cooling for either counter or 
parallel flow. Prints the coolant mass « 
flow rate and pressure. 


Calculates the specific heat at constant 
pressure as a function of temperature for 
routine HEAT and the specific heat ratio as 
a function of Mach number of routines 
RIMAIN and HEAT. 


Computes Kroneker Delta for matrix solutions. 


Determines the real roots of a third order 
equation. 


Computes the matrix inversion and accompanying 


solution for linear equations. 


Computes the distance between two given 


+ tac 
coordinates. 


Selects the radiation view factor data from 


the input data and forms the radiation matrix. 


Calculates a pressure balance and surface 
temperatures for a film cooled surface. 
Prints the coolant mass flow rate and 
pressure, 
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A-122 


A-40 


A-139 


A-128 


A-196 
A-178 
A-199 
A-215 


A-202 


A-134 
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ROUTINE 
NAGE 


FINDIT 


FLOWNO 


GATHER 


HEAT 


HEATTR. 


HOTPT 


MATPOS 


MATSYM 


MINV 


ORDER. 


PICMAX 


PICMIN 


PLANCK 


PLOTL 


TABLE I (Cont) 


RASIC FUNCTION 


Locates the axinl positions of two nodes with 
respect to the axial positions of the confipura- 
tions surfaces. 


Checks the resulting pressure for surface 
coolince against the stotic pressure wions 
the cooled suriace. Prints a warnir. 
diagnas vic foe greas of nu cooling fier. 
Combines data “~om routines AXAREA and 
COPLOW to form ce compressible flow 
pronerties for cas consiguration flow 


STTOaMS , 


Calculates the initial conditions necessary 
for the boundary iuyer calculations. 


Calculates the convection heat transfer 
coefficiont. 


Routine to control the calculation of 
surface boundary izyers and heut transfer 
coefficients. 


Checks to ensur: thet all the elements for 
convection and radiation matrices are 
positive. 


Completes the matrices for the convection 
heat paths and the radiction view factor 
area product, 


Inverts the "Script F" matrix, 
Arranges two-dimensional coordinate arrays 
into increasing order according to a 


selected variable. 


Determines the maximum value of an array of 
up to twenty numbers. 


Determines the minimum value of an array of 
up to twenty numbers. 


PLANCK Black uvoedy function. 
Generates a plot of radiant intensity vs. 


wavelength for both the non attenuated and 
the atmospheric attenvated radiation data. 


PAGE 
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A-1I¢ 


A127 


A-153 


A-159 


A-149 


A-197 


A-198 


A-222 


A-190 


A-191 


A-219 


A-~206 
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ROUTINE 
NAME 


PRND7T1 


PRNDTL 


PSOLN 


QUAD 


REDI 


REDII 


REMAIN 


SETFLO 


SHADOW 


TABLE, 1 (Cont) 


BASIC FUNCTION 


Calculate the Prandtl number as a function of 
termperature for rgutine CONFLM. 


j 
Calculates the Pradtl number as a function of 
temperature for routines TURBLT and RIMAIN. 


Generates surface cooling flow supply and 
discharge characteristics and solves these 
simultaneously for calculation of the available 
surface cooling wass flow rate and pressure. 


Determines the roots of a second order 
equation, 


Computes the spectral radiant intensity 
(watts/Steradian) for each engine off-axis 
angle. 


Computes and prints geometric view factors 
for a detector to the configurations 
surface nodes as. a function of engine off- 
axis angle. 


Control routine for the fluid flow, heat 
transfer, surface cooling and surface force 
factor calculations. Prints the fluid flow 
information. 


Determines the intersection of a shadowing 
body with the "cone of vision" that a point 
on one node secs another node. 


Used to print the heading informiution for 
routine REDII. 


Prints the surface radiation information. 


Computes internal radiation interchange 
factors for the thermal analyzer. 


Determines the pressure balance for the 
surface cooling routines if a cooling mass 
flow rate is specified. 


Describes the configuration surfaces uscd for 
shadowing in routine VIEW in a manner compatible 
for shadowing surface description of routine 
REDIT, 


A-120 


A-177 


A-220 


A-210 


A-110 


A-189 


A-216 


A-224 


A-204 


A-133, A-138 


A-209 


ROUTINE 
NAME 


SIGNIR 


STORE 


“STOREY 


STRATE 


SUMAT1 


TARLEL 


TAIL 


TAIL1L 


TAILLE 


TAS1R 


TESTI 


TEST2 


TESTN 


TABLE 1 (Cont) 


BASIS. FUNCT. ON 
The programs main routine. Controls the 
flow throughout the entire noresram, 


Temporarily stores information necessary 
for everage heat *rensser coctTicient 


5 


ep leulat ons. 


Temporarily stores information vecessarv 
for force factor caiculations. 


Computes solution utilizing Linear inter- 
polation routine used in ceseribing the 
configuracion surface ata Cor boundary 
layer calculations. 


Summates and prints tne adiant intensity 
for selected wavelennsch bande, 


Computes solution utilizing second order 
interpolation of dita and provides the 
first derivetive of the intermolated 
function, 


Interface parameter generation and 1/0 
program control routanc. 


Metal signature calculation control routing. 


I/O and interface parameter calculation 
control routine. 


Controls program flow snd data involved in 
calculating the radiation interchange factor 


Matrix, conduction matrix, performing the 
thermal balance and printins the wall 
temperature for cach ade, 


Tests the possibility of one node being 


shadowed by a second node which has the same 


axial location. 


Uetermines whether » point lies between two 


other points. 


Computes the last computed value of geometric 


view factor with the previously computed 


value to determine whether the view factor 


Needs to be reci:lculated using finer 


divisions of the nodes. This routine also 


Sets these divisions. 


A-125 


A-124 


A-121 


A-225 


A-205 


A-208 


A-218 


A-193 


A-180 


A-181 


A-192 


TABLE 1 (Cont) 


: ROUTINE 
NAME BASIC FUNCTION PAGE 
TITLA Controls output printing of the title page. A-4] 
TRANCL Calculates a pressure balance and wall temp- 
eratures for a transpiration cooled surface. 
Prints the coolant mass flow rate and 
pressure. A-130 
TURBLT Calculates and prints the boundary layer and 
heat transfer information as a function of 
configt:ration axial distance. A-154 
VECANG Calculutes the argple of a node relative to 
the configuration center line. A-217 
VIEW Calculates and print: for each configuration 
node, internal geomet.ic veiw factors. A-162 
WRIT14 Transfers a two-dimensional array into a 
“ single dimensional array. A~201 
4? ic 
\ 
{ 
; 4 
t . 
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* { 
1 23 “ 


oo The equivalent Black Body Area and Temperature are extracted from the 
eo data generated by SIGIR. The first step is to determine the 
eae equivalent Black Body Temperature. This is accomplished by searching 
the generated spectral data for the maximun intensity, figure 4, and 
applving the Wien Displacement Law: 


The equivalent Black Body Area is found by dividing the integrated 
signature obtained from SIGNIR by the Radiant Exitance over the band cf 
intere'st. 


iE 


te 42 ie 
Agg = { Ma(STGNIR) / [ Me(Tpp) 
Ay se 


Figure 5 shows a typical correlation hetween the predicted spectral 
signature and the sncctral signature of the equivilent hlack body. 

These resultant temperatures and arcas are the only neramceters transferred 
hack to ASDIR. 
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PROGRAM: 


FUNCTION: 


INPUT: 


OUTPUT: 


SUBROUTINES : 


DESCRIPTION: 


PLMDEF (Page A-81) 


Plume definition .and gas data array development. 


Engine nozzle exit flow properties (velocity, pressure, 
temperature, and species [H20, C02] concentrations) as 
‘a function of radius. 


Plume gas data array. At this point, velocity (ft/sec) 
and pressure (atm x 107%) are multiplexed into common 
registers. 


FLINP 


The gas dynamic calculations made for the plume in this 
program, figure 6, are applicable to axisymmetric 
turbojet and turbofan engine with mixed or separate 
co-annular co-plonar exhaust jets. The finite difference 
calculations were developed to express the conservation 
(mass, momentim, energy, and species) laws for numerical 
integration by General Electric in their early SCORPIO 
work reported in reference [5]. Two effective kinematic 
viscosities were retained, one for physical mixing and - 
the other for thermal mixing. The conservation laws 
were simplified for axisymmetric flow whereby all 
gradients in rotation about the centerline were not to 
zero, and all radial components of momentun were 
neglected as being small. Second axial derivatives and 
axial gradients divided by the radius were also neglected 
in momentum. Thermal radiative heat transfer and thermal 
generation were also neglected. In the derivation of the 
energy conservation and transfer functions the effects of 
the second viscosity, the square of axial gradients, 
dependence upon radial velocity, and gradients divided by 
the radius were minimized. The radial coordinate system 
was transformed into a stream function cucrdinate system 
after Von Mises. The transformation functions end the 
transformed conservation relations to which the finite 
difference procedure was applied are as follows: 


Stream Function 


wherein: 2¢ « eel oo = 


or r ” @x 
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4,0 "PLMDEF | 


} INITIALIZE _ 
| 7 PRESET 


“T CALC: RADIT FOR ‘STREAMLINES 
LOAD OUTPUT ARRAY FCR STATION 
eee as enseek eee eens Sennen ene ene 
LF FINAL STATTON E Decca 


<A 


i 


RESET FOR NEW STATION 


SEARCH FOR SHEAR .LAYERS 
"ADJUST | NCREMENTAT ION 
“ADD AMBIENT RADI] AS NEEDED, ; 


RESET MIXER FOR FIRST SHEAR LAYER 


ANCE MIXER FOR SECOND SHEAR LAYER. 


F STABILITY TESY SHOWS UNSTABLE 


LOAD PLUME GAS: DAT'A ARRAY 


F STATION STAGING CRITERIA NOT SATISIFIED 


ne fe 


TO RADII: CALCULATION 


4 
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Figure 6. PLMDEF Flow Diagram 7 d 


CONTINUITY: 
for the stream function: 


ov ov  vo'w ay. ay 
=n Or ox” Orox ~ ox or” 


and after the transformation: 


nS je on _9/ pyr =O 
ox{ pur d¥ | pur 


ov, ¥ av pudx 
ENERGY 
aT, [pul \(_a}\ aul? , 2 (o%ue? ar’ d 
ar (ours (34) (3 + al a¥) * pe ax 1 
4 
ur? 1 yi 
+ ey) pan og 1 
|: 
DIFFUSION: 4 
i! d 
oA 
bare = Bag peur? 284) 4 
ox = Yo¥, Y * ay} 
where: | 
T = gas static temperature or thermal subscript 
u = gas axial velocity oe 
\ 


v = gas radial velocity 
Prandtl number ~« c pes =] 


Schmidt number = €/D = 1; in which D is the binary ‘| 
diffusion coefficient 


€ = effective kinematic viscogiiy | 


oO, ™ species concentration 


ac 
== is an assumed radial dependence for the gradient 


of the specific heat of a mixture of gases = c/t. 


Other quantities and terms are conventional. 


The effective viscosity adopted for ASDIR includes the Donaldson and 
Gray treatment for the compressible frce mixing of two dissimilar gases 
as reported in the ALAA journal, Val 4, no, 11, in November 1906 and as 
applied by General Electric in SCORFIO-11 [5]. 


Overall, the plume model appears to reflec’ the cssential features of 
measured plume data in those rewions where the gas is hot. The regions 
where the model tends to lose accuracy are toward the outer edges where 
‘he gas shear appears to be too small thus allowing the plume to spread 
‘xcessively. However, the plume temperatures and species concentration 


are essentially ambient in wmiue so that the radiance und transmittance 
errors are small. 


Viscosity Model 


The choice of the turbulent viscosity, .¢, is critical to the 
P plume flow field. Past models, available in the iiterature (Schlicting, 
‘ Abramovich, Varren and Applied Science Laboratories), show the turbulent qt, 
“viscosity as a function of axial distance and velocitics in the mixing ; 
"region. Preliminary results for the flow field showed that the analysis ! 
using these models did. not agree with jet engine test data. A modified | 
temperature dependent model for viscosity was formulated to include the | 
temperature effects of hot flow. 


The model being used in this analysis is: 


c= xe{Uh = ULI LENS: (T,/T2) 7] +y } 
’ ts 
where 
y=13 M,< 0.6 } 
74 
y 1.6/1 +M) 3 M > 0.6 
4 
The value for y is from a paper by Donaldson and Gray. The ; 
- r value of CNST depends on the nozzle exit Mach number (M,) and the centerline a 
. & ‘mperature. (Sec Page A-85). 
Stability and Print Criteria 


In order to insure stability, the axial step size, Ax, must be 
fe kept small. It was found that the system became unstable when the 
" influence parameter 


1 
pee 


hs ez = {Axtceurr?=p?}/{2-y?+T*-r7} > 0.20 | ; 


“ 
aa) 
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To achieve a solution in a reasonable length of time, the 
step size Ax should be as large as possible. Accordingly, czm, the 
: maximum value of cz for each axial location is examined. If czm is 
too great, Ax is reduced for ensuing calculations. If czm is too 
small, Ax is increased for ensuing calculations. 


Flow parameters are not entered into the data array at every 
value of Ax calculated. The calculation step sizes are too small. 
Tabulation occurs whenever the number of Ax increments from the previous 
calculation exceeds NBS or the sum of the Ax's exceed Axp. Axp is 
incremented after each tabulation such that 


AXPN = XPg (0.5N) 


Excessive static pressure at the nozzle exit plune is treated in a 
pseudo convergent-divergent nozzle manner which reduces the pressure to 
ambient value while accelerating and expanding the plume core to 
Supersonic velocity over an axial distance of one nozzle diameter. 


The program develops a plume gas data array containing the velocity, 
pressure, temperature, COz2 concentration, H20 cencentration, and radius 
for up to 30 radii over an axial span of 49 stations. The velocity and 


pressure are multiplexed to share cammon registers. (See PLMDM for 
demultiplexing) . 
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PROGRAM: 


FUNCTION: 


INPUT: 


OUTPUT: 


SUBROUTINES: 


DESCRIPTION: 


FLINP (Page A-87) 


Previde a vehicle by which the required calculation 
pavaneters, nozzle exit flow quantities und data 
initialization can be read into ASDTR. 


~All control and computational variables in the PLUMIN 


name List. 


Computationally compatible cngine nozzle exit plane 
Flow quantity data. Nozzle and wumpient data are printed. 


CHEM, THRUST 


The flow diagram of Program FLINP is shown in Figure 7. 
The function of this program is to initialize ensrine 
operational input and control quantities. The control 
function is derived from the input data. If the engine 
operational input quantities are not of the nroper forn, 
certain operntions and additional reads are invoked to 
put the input into a fona compatible with the PLMDEF 
program. Details can be found in ASDIR-II, Volune I, 
Users Manual. 


In general, if PA is contained in namelist PLUMIN, the 
Mil Std 210 standard atmosphere in subroutine THRUST is 
bypassed, The inference is that an off-standard or test 
facility ambient condition is being simulated and the 
other ambicnt quantities are required in PLUMIN. 
Similarly, this program will bypass the default engine 
sample case if a nozzle exit velocity U8(1) is piven 
greater than 1 ft per sec or if an ambient velocity UA 
is given greater than 11 ft per sec. Further, if U8(1) 
i, greater than L the volume for T8T, P8, PY, XCO2Z, XH20 
will be accepted as piven in PLUMIN. Again, if xCU2 has 
been omitted or given a small value (4.c. « XCO2ZA), the 
XCO2, XH26, R, GAMMA, etc. will be calculnted in CTEM 
for the fuel type (TANE) and stoichiometric cquivalent 
ratio (TQR). Finally, if a gas quantity distribution for 
U8, T8T, XCO2, and XH20 is piven in PLUMIN, they will be 
accepted unchanged since the load urray key is taken from 
a test on U8&(2). 
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Figure 7. FLINP Flow Diagram 
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PROGRAM: 


: FUNCTION: 


TNPUT: 


OUTPUT: 


SUBROUTINES: 


DESCRIPTION: 


CHEM (Page A-91) 


Generate the post combusticn gas description (gas 
constant, specific heat and speci “ic heat ratio, 
molecular weight, species conserv ition, etc) for a given 
fuel type, TANT, stoichiometric equivalence ratio, TQR. 


TANE, FOR 


Fuel to air ratio, gas constant, specific heat, gamme, 
CO2 concentration, and 20 concentration. 


None 


This program calculates the gas properties for the 
products of combustion with ideal combustion. Since 
combustion products properties vary only slightly with 
pressure and tend to chemically freeze neur 2500°F in 
typical nozzle flow, the precision of this simple college 
chemistry «approach is adequate and the fist calculation 
makes it a worthwhile approach. 
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PROGRAM: 


FUNCTION: 


INPUT: 


OUTPUT: 


SUBROUTINES : 


DESCRIPTION: 


THRUST (Paye A-93) 


Provide a vehicle by which the required calculation 
parameters and engine operating quantities can be read 
into ASDIR. 


All control and computational variables in the POWER 
namelist. 


Comp.cationally compatible engine nozzle exit plane 
fl-w quantities not provided in FLINP input namelist 
PLUMIN. Flight condition data is printed. 


CHEM 


The flow diagram of subroutine THRUST is shown in figure 
8. The function of this program is to read atmospheric, 
engine type, data type code, and engine specific 
performance data. The program control functions described 
for Program FLINP apply throughout this program. 
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| IF (u8(1) > a 
i” IF (UA > 11) 
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ee ae 
CALC RAM 
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| LOAD DATA ARRAY 


PRINT FLIGHT CONDITTONS | 


Figure 8. THRUST Flow Diagram 
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PROGRAM: 


FUNCTION: 


INPUT: 


OUTPUT: 


SUBROUTINES: 


DESCRIPTION: 


PLMOM (Page A-98) 


To manage the plume gas data. An input control 
parameter, KDATA, 


gas data. 


Plume gas data and the control code KDATA. 


Plume gas data output to cards, tape, printer, or 
printer plot in any combination. 


PLMPLT, PLMPRNT 


Each digit of the five digit integer KDATA cmtrols an 
I/O function. KDATA is represented by the valuc: 


KDATA = ABCDE 
Tf A: 0, a new plume will be calculated. If A= 1%, 


an old plume previously filed will be brought in and 
plume calculations will be bypassed. If A> 1, plume 


gas data will be 


If B = 1*, the plume gas duta array is filed in the 
computer for later or repeated analysis. If B #1, the 


tiling operation 


If C = 1%, the plume gas data array is punched into 


output cards for 


the punch operation is bypassed. 


If D = 1, the plume gas data array is printed in output 
by subroutine PLMPRNT. If D # 1, the print operation is , 


bypassed. 


The subroutine PIMPLT is called with the final valuc E 
in the KDATA code. The i value has a different pattern 
than exhibited by A through D, see PLYLT. 


controls the dispatch of the plume 


read from input cards. 


is bypassed. 


later or repeated analysis. If C #1, 


Tn addition to breaking down the code carricd by KDATA, 
this program auds a layer of ambicnt data all around 


the plume data. 


The ambient data blanket is necessary 


to contain the interpolation exercises performed in the 
ray tracing calculations. The program then searches the 
plume gas data array to locate the tip of the core and 
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record such other parameters such as plume length AL, plume diameter 
REND, plume edge tanzcnt TNR, etc. The minimum output RB, XC, REND, 
and AL is printed. 


The PLMDM flow chart is shown in Figure 9. 


*This code represents a machine dependent opcration and may require 
fur ther deveionment for satisfactory oreration. 
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Figure 9. PLMOM Flow Diagram 
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OUTPUT: 


SUBROUTINE : 


DESCRIPTION: 


PLMPRNT (See Page A-100) 


Print the contents of the plume gas data array on 2 
the output printer. 


The plume gas data array 
The plume gas data 
None 


Print the plume gas data array. The plume gas data is 
printed in the following sequence: 


- general caption 

- plume station number and location 

- centerline and edge Mach numbers 

- radius, velocity, pressure, temperature, and C02 
and H20 concentration 


Units are printed in the captions. The quantities printed 
are all self explanatory except the edge Mach number. 

The edge Mach number is that valuc at the outer edge of 
the innermost shear plane. This quantity will tend to 
oscilate due to the stepping and clipping process which 
occurs in the finite differencing procedure. 


eb nie ee chia A 


Dee ent ara a Mice anh ata ett ante 


PROGRAM: 


FUNCTION: 


INPUT: 


OUTPUT: 


SUBROUTINES : 


DESCRIPTION: 


PLMPLT (See Page A-101)- 


Plot selected plume gas data quantities on the 
output printer. 


The plume gas data array and the final digit of the 
KDATA control code. 


Selected arguments of the plume gas data plotted on 
the output printer. 


None 


The output printer has 136 characters horizontally and 
more than 60 lines vertically. The grid represented 

by lines and characters are assigned one of the letters 
in the word plume (P L UME) with the additio of # 

to permit the plotting of six ranges. A space is 
assigned the first character portion after a change in 
range, The nozzle exit plane value establishes the 
separator between # for values greater than the nozzle 
exit value and P for values between 80% and 100% of the 
nozzle exit value. Values in the 60% to 80% range are 
assigned L. Values between 40% and 60% are assigned 

U, etc. Each line represents a radius, and the 
interstation distance accumulates a aumber of characters. 
With two way interpolation a printer line is compiled 
and printed. After the data has been plotted, the axial 
scale is compiled and printed. 


The code designating the plume gas data to be plotted is 
provided by the single digit value remaining in KDATA 
after having been stripped in program PLMDM. Up to 

four plots can be generated. In their plotting sequence, 
they are temperature, X02 concentration, H20 concentration, 
and plume velocity. The appropriate captions are 
included. All four data arguments are plotted if the E 
value of KDATA is given as 0. If E = 1, no plots are 
printed. If E = 2, the last three are printed, if E = 3, 
the last two are printed, if E = 4, the last plot only 
is printed, and if E is greater than 4, the program will 
fail. For the general breakdown of KDATA, see program 
PLMDM. 


The final activity accomplished in program PLMPLT is the 
de-multiplexing of the velocity and pressure. The 
pressure is brought up from the back of each velocity- 
pressure register and replaces the velocity. The 


velocity is no longer retained in the plume gas data 
array. 
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PROGRAM: 


FUNCTION: 


INFUT: 


OUTPUT: 


SUBROUTINES: 


DESCRIPT TON: 


ALPLUM (See Page A-47) 


Calculate total apparent spectral radiance of 
exhaust plume and hot cngine parts us reccived by 

an observer with plume and atirospheric cmissions and 
attenuation, including background contributions 
considerec. 


Flight conditions, engine hot parts and plume 
parumeters. 


Engine/Plume IR Signature 
RAYCAL, PLURAY, PLUSIG 


The flow diagram of program ALPLUM is shown in figure 
10. The prime Function of this routine is to be the 
controlling program to allow calculation of aircraft 
infrared signatures. Description of the working 
programs controlled by ALPLUM are contained under the 
individual program titles. 
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PROGRAM: RAYCAL (Page A-50) 


FUNCTION: To segment plume into discrete areas for signature ‘ 


analysis in other routines. 


INPUT: Plume ard engine geometry. 

OUTPUT : Array of plume scement parameters. 
SUBROUTINES: None 

DESCRIFTION: 


RAYCAL, being basically geometric, is a relutively 
simple program. Nowever, since signature prediction 
uccuracy is deyenden! upon the geometric representation 
of the engine nozzle and the plume, graphical detail 
will be presented for this routine. 


The engine nozzle and plume are characterized as a series he 
of planar area segments each representing the termination 

of a ray or beam of radiated energy. Lach segment 
represents cither free (background) space or a black 
surface. The cross-s¢tional shape and area of each ray 
is defined by the seguentation scheme and is constant 
along the ray centeriine. The segmentation scheme takes 
full advantage of symmetry us described with reference 
to the six accompanyiny figures. 
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Rectangular coordinates provide, figure 11, the frame- 
work for the ray definition. The X axis represents the 
engine nozzle centerline. The D axis, in the nozzle 
exit and perpendiculur to the centerline X, is free to 

' rotate about X so that the plane D-X will contain the 
i centerline ray. The lateral axis, Z, is normal to both 
" D and X. The angle of view, known us the aspect angle 
or ASPDEG, A8, is meusured from the centerline X and 
: lies in the D-X plane as shown in figure 12. 


the ray centerline on plane D-Z from the Z axis. 


Dimensions on D represent the normal intercept of uy 
4 
Dimensions on Z represent the lateral intercept of re 
the ray centerline on plane N-Z from the ) axis, oy 
Dimensions on X represent the axial intercept of the 
ray centerline on plane X-Z from the Z axis. 


re 

The simplified plume geometry adopted for use in this ; 4 

. program is shown in figure 13. The plume is analyzed in ihe 
= three independent but adjoining regions: | 


(1) The plume core is segmented in the X-Z plane into a 
number of radial and angular segments shown in figure 14, 
“ ae Segmentation begins at the extremity of the nozzle exit 


projection on the X-Z plane and continues to the extremity of 
| the core. 


(2) The plume skirt is segmented in the X-Z plane into a 
nunber of radial and angular segments shown in figure 15. 


Segmentation begins at the extremity of the core and continues 
to the extremity of the plume. 


(3) The engine nozzle is segmented in the D-Z plane into a 
number of radial and angular segments shown in figure 16, 


Each planar segment so gencrated also defines a ray segment. The D-2Z 
or D-X ray segment centerline intercept is centrally located, see 
asterisks, in each planar segment. Further, each segment is considered 
to be the ray segments own projection on the D-Z and D-X planes from 
which the ray segment cross section area, RAR, is calculated. 


The ray segments are generated on the first call to RAYCAL. Successive 
calls simply transfer previously generated data. 
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FUNCTION: 


INPUT: 


OUTPUT: 


SUBROUTINES : 


DESCRIPTION: 


PLURAY (Page A-54) 


Generate the parameters describing the ray which 
are required in the signature analysis routines. 


Ray starting geometry, plume gas parameters. 


Ray geometry, gas partial pressures and gas static 
temperature arrays. 


INTERP, START 


By referencing the plume gas data and the ray 
definition parameters, this routine fills the pressure, 
temperature, and length urray with the proper data for 
ray Signature analysis, see figure 17. 
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PROGRAM: 


FUNCTION: 


INPUT: 


OUTPUT: 


SUBROUTINES: 


DESCRIPTION: 


START (Page A-57) 


To calculate the starting and ending points of 
a ray through the plume. 


Plume geometry, ray intersections 
Starting and ending dimensions of ray 
None 


The starting and ending points are banded by 
engine geometry, engine plug geometry and flow 
field boundaries. 
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OUTPUT: 


SUBROUTINES : 


DESCRIPTION: 


INTERP (Page A-61) 


Interpolate on gas parameters to find values at a 
given point y ‘lin the plume. 


Gas parameters and coordinates of point in plume. 
Concentration and temperature at input point. 
Nor, 


This program interpolates between stations along the 
centerline of the plume and along the radial dimension 
of the plume to find the static temperature, mole 
concentrations of water vapor and carbon dioxide and 
gas pressure at the input coordinate. 


PROGRAM: PLUSIG (Page A-63) 


FUNCTION: Evaluate the spectral signature and transmittance 
of a nonuniform gaseous path of NT segments. 


INPUT: Ray properties (temp, pressures and lengths), and 
frequency interval. 


OUTPUT: Spectral signature of gaseous ray and spectral 
transmittance of ray. 


SUBROUTINES: ATMOS, KBCAL, PLANCK, INTERPO, SETTAU, TAUCAL, 
TAUN20, ERF, DAH20, DCO2 


DESCRIPTION: The flow diagram for PLUSIG is shown in figure 18. 
The prime function of this program is to calculate 
the spectral signutures of nonuniform gas paths with 
various end conditions. The gaseous ray is divided 
into NI nonuniform segments with each segment defined 
by the partial pressures of H70 (PHj), CO» (PC;) and 
broadening pas (PN;), length tD;) and static | 
temperature ([T;). k brief discussion of spectrai 
radiators follows: 


| 
Continuum radiation is obtained from liquids and 
solids because the molecules of the substance are in 
direct contact and exert forces on each other causing 
all levels of molecular emissions to occur. However, J 
in the case of gas molecules, ihe intermolecular i 
spacing is large and the free molecular vibrations 
and rotations result in discrete spectral emissions. 


The intensity of these discrete gaseous emissiors 
obeys the PLANCK radiation law as modified by Kirchoff's 
law. However, the emissivity of the gas is spectrally 
dependent and radiation is emitted or absorbed only in 
narrow wavelength intervals called spectral lines. 

These spectral lines are clustered into tightly spaced 
sets called bands. The emissivity of the gas is zero 4 
outside these bands and dependent upon line spacing e 3 
and intensity within the bands. 
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The inhomogeneous gas band-medel expressions used a 
in this program were developed at the University of ae 
t Michigan and are based on the Curtis-Godson a 
approximation. Rather than attempting to fully explain 
the excellent work perfomed at the University of 
Mic''igan Infrareu and Optics Laboratory (now ERIM) in 


ee 


-. the uveas of gaseous emission\and absorptians, see 4 
oc references [5, 6, 7, 8, 9, 10]: ae 
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A simplified discussion of the function of the program PLUSIG 
follows: <A ray of NI segments, NA of which are atmospheric, of 
cross sectional area A is entered into the program, figure 19. 

Each of the NI-NA segments, representing the exhaust plume gas, are 
described by the gaseous partial pressures, temperatures, and 
lengths. The atmospheric segments are front loaded in the ray array 
through the use of the subroutine ATMOS and the range between, and 
altitudes of the observer and vehicle. The transmissions of each 
Yay segment are calculated in PLUSIG using subroutines KBCAL, 
INTLERPO, SETTAU, TAUCAL, TAUN20, ERF and data contained in Data 
Blocks DAH20 and DCO02. 


Figure 19 depicts the information available for a given 
frequency interval after the exercise of the transmittance 
calculation subroutines. The signature of the ray is found from 


NT 
Fe eae Se @ 


Jra. 
ray i=] ©Y 1 vi 


ev 
where Ley is the PLANCK function and a ei is the absorptance when 


i 
a; = (1 - Ty) ol t 


v1 jel vj 


The total transmittance, t, for the gaseous puth containing N1' 
elements is extractable, 


NT 
il ts 
ray i=l i 


Thus the spectral transmittance can be used to determine the 
piume and atmospheric attenuation of a black surface at the 
temminus of the ray, figure 20. 
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Figure 20. Ray Geometry with Temninating Black Surface 
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PROGRAM: ATMOS (Page A-69) 


FUNCTION: Generate temperatures and partial pressures along 
: an atmospheric path of NA elements and front load 
= the ray pressure and temperature array with the 
atmospheric data. 


INPUT: Range, altitude of plume and observer, aanbee of 
atmospheric segments 

OUTPUT: Hz0, CO2 and N2 Pressures 

SUBROUI!'INES: None 

DESCRIPTION: Given the geometry of figure 21, the altitude at 


point r is generated and the temperature and 
pressures are found from the empirical curve fits 
used in the computer code. References [8, 11, 12, 
13]. 
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PROGRAM: 


FUNCTION: 


INPUT: 


OUTPUT: 


SUBROUTINE: 


DESCRIPTICN: 


or 


SETTAU (Page A-70) 


Insert pressure broadening effect of broadening 
gas by modifying effective gas pressures. 


H20 or CO2 and N2 gas pressures 
Modified gas pressures (H20-N2 and CO2-N2 mixtures) 
None 


References [6, 7, 8, 9, 10] 
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PROGRAM: TAUCA. (Page A-71) 


FUNCTION: Calculate the gaseous transmittance in a given ° 


frequency interval. 


INPUT: Pressuves, temperatures, gs parameters (¥ & 8) 
QUTPUT : Gascous Transmittances, TCV) 
SUBROUTINES: INTERPO 
DESCRIPTION: Calculate transmittance of gascous substance over a 
given frequency interval using overlapping line 
upproximat tons. 
é 
t(v) = exp[-Ry°f] 4 
| 
Where Pe is the effective value of and f is the 
ie. Ladenbure-Reiche function. References [6, 10]. te 
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PROGRAM: 


FUNCTION: 


INPUT: 


QUTPU?: 


SUBROUTINES : 


DESCRIPTION: 


KBCAL (Page A-73) 


To look-up and interpolate on the band model 
parameters for a given frequency. 


Frequency, gas type 
K, 8, temp and number of temperatures 


None 


The average line width to spacing ratio multiplied 
by 1/2 1 called 8, and the average line strength to 
spacing ratio called K are generated in this 
subroutine. These parameters are found at a given 
input frequency for H70 and C02 over th temperati ‘2 
range from 300°K to 3000°K. 
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_ 


OUTPUT : 


SUBROUTINES : 


DESCRIPTION: 


Finction PLANCK (Pase A-75) 


Calculate PLANCK's blackbody function 
Frequency (em7!), temperature (°K). 
Watts/Steradian/on> 

None 


From the classical expression for the radiant 
emittance into a Lanbertian hemisphere 


= 2u(KT)* x" 
Mey ale X)- | 


where 
x = hv/eT 


the spectral radiant sterance for a given frequency 
and temperature can be found 


les 2he?v3/ (exp (x) -1) 


See Reference [14]. 


64 


he 


went ee ee i 


cabal. ad see tie wat ah Oh 


— il 


t 
{ 

‘ 
4 
i 
4 


PROGRAM: 


FUNCTION: 


INPUT: 


OUTPUT: 


SUBROUTINES: 


DESCRIPTION: 


Function INTERPO (Page A-76) 


Interpolate over the N dimensional set (X, Y) to 
find the value at a given point (ARG). 


X, Y, N, ARG 
INTERPO 
None 


See Puge A-76, 
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cake PROGRAM! TAUN20 (Parc Ae77) 


FUNCTION: Calculate transm.ttunce of NoQ in the atrosphere : 
over colecte! froaucncy interwd:, 


INPUT 


Frequency, atmeswheric Vat. 


OUITUT: Transmittance of N90 in atmosnacre. 


SUBROUP INES: RF, INTGRFO 


: 
3 
a 4, ‘| 
NIGCRIPTION: See Page A-77. : 
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PROGRAM: Function ERF (Page A-78) 


FUNCTICN: To evaluate the error function, erf(x). 
INPUT: X 
OUTPUT: ERF (x) 
SUBROUTINES; None 
&. DESCRIPTICN: Evaluates the error function. 
x 


Serf(x) = 2//r fempc-t2,ae 


12. 


13. 
14, 


ee a ee er a ere nad 


Re Teronces 


Ludwig, C. B., ‘Measurements of the Curves-::°-Growth of Hot Water Vapor", 
APPLIED OPTICS, Vol. 10, No. 5, Mav 1971, nos 1987-1975. 


Huffaker, R. M. and M. J. Dask, "A General Program for the Calculation 
of Radiation from an Inhomogeneous, Non-isoxciis, Nen-isethermal Rocket 
ixhaust Plume". Aero-Astrodynamics Leberetory, Cenres C. Vershall Space 
Flight Center, "National Agronauties and “rece Administration, NASA TM 
N-S5622, dune 19, 1967. 


Parcells, R. F., A. N. Blaser and T. ™. Mouricsen, “A rerart Engine 
Infrared Signature Prediction Program (GUGNTR-2). USTR'S MANUAG', Vought 
\eronautics, LTV Aerospace Cerporcion, IM Ne, 2-Ss8ReO/SR-2997, 

25 June 1975. 


Parcelis, R. F., A. N. Blaser and T. T. Mouritsen, "Aircraft Engine 
{Infrared Si ignature Prediction Program 'STGSIR}, USER'S MANUAL, Vought 
Aeronautics, iTV Acrospice Corposution, TR Nor R2-55300/1R-2947, 

30 Sept 1971. 


Wilton, M. E., ''An Analytical Technique for Spectral Calculatior of 
Radiation in a Plume including Ghstructions (SCORPTQ)", General Electric 
Report No. R7QAFG299, April 1970. 


Lindquist, G., "How tc Use FS + TAUS'', University of Michigan, MEMO, 
28 Jan 1970 (see ASDIRK-JI, Volume TIT, Reference Documentation). 


Lindquist, G., "Additions and Changes to US", University of Michigan, MEMO, 
17 April 1970 (se. Volume TIf). 


Lindquist, G., "How to Use the Program Plumes", University of Michigan, 
MEM), 22 July 1°70 (see Volume TIL). 


Simmons, F. S., and G. Lindquist, “IR Radiation from Recket Exhaust Plumes 
in Goace, University of Michipan, Ri2i2-1-", Oct 1970 (see Volume TII). 


Lindquist, G., and F. S. Simmons, "A Band Model Formulation for Very 
Nonuniform PATAS"”, J. Quant. Spectral Radiant Transfer, Vol 12, Pps 807- 
820, Pergamon Press 1972, 


A Reema cecins one 5 -aTTNm gan 15 
d Transmission anc Backs? ound Radiation by 
x q 
t 


nt No. GLEDI99, Dec 1961. “(n401923 


tno} 


U.S. Standard Atmosphere Supplements, 1966, ESSA, NASA, USAF, 1966. 
Byers, H. B., General Meteorology, McCraw-Hall Book Co., 1959. 


Handbook of Military Infrared Technology, Office of Naval Research, Dept. 
oF the Navy > LO5. 


68 


pedi A EE ee ue 


mz 


APPENDIX 


ASDIR-II 


PROGRAM LISTING 


St mens, cate = 


APPENDIX 


This appendix contains a listing of the Aeronautical Systems Division's 
Infrared Signature Predicting Model (ASDIR). The computer code was 
developed in the structural form shown in figure 22. The overlays are 
organized in three major functional groupings with a fourth overlay 


dedicated to data initialization. The prime fumctions of each overlay 
are: 


Overlay 

(0,0) Program Control 

(1,0) Data Input § Initializatioa 

(2,0) Plume Signature Calculations 

(3,0) Data Initialization 

(4,0) Plume Gas Dynamics 

(5,0) Plume Data 1/0 

(6,0) Spectral Filter Data Initialization 
(7,0) Spatial Output 

(11,0) Hot Parts Temperature Calculation 
(12,0) View Factor Calculations 

(13,0) Hot Part Radiation Calculations 
(14,0) hot Part Data Conirol as 
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